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We explore scenarios where the highest energy cosmic rays (HECR) are produced by new particle 
physics near the grand unification scale. Using detailed numerical simulations of extragalactic 
cosmic and 7-ray propagation, we show the existence of a significant parameter space for which 
such scenarios are consistent with all observational constraints. An average fraction of ~ 10% 
7-rays in the total cosmic ray flux around lOEeV (10 19 eV) would imply both a non-acceleration 
origin of HECR and a large scale extragalactic magnetic field < 10~ n G. Proposed observatories for 
ultra-high energy cosmic rays should be able to test for this signature. 
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The HECR events observed above 100 EeV are 
difficult to explain within conventional models involv- 
ing first order Fermi acceleration of charged particles at 
astrophysical shocks pJ. It is hard to accelerate pro- 
tons and heavy nuclei up to such energies even in the 
most powerful astrophysical objects 0], like radio galax- 
ies and active galactic nuclei. Also, nucleons above 
~ 70 EeV undergo photopion production on the cosmic 
microwave background (CMB), which is known as the 
Greisen-Zatsepin-Kuzmin (GZK) effect B and limits the 
distance to possible sources to less than ~ 100 Mpc ||. 
Heavy nuclei are photodisintegrated in the CMB within 
a few Mpc @] . There are no obvious astronomical sources 
within ~ 100 Mpc of the Earth. 

A way around these difficulties is to suppose the HECR 
are created directly at energies comparable to or exceed- 
ing the observed ones rather than being accelerated from 
lower energies. In the current versions of such "top- 
down" (TD) scenarios, predominantly 7-rays and neutri- 
nos are initially produced at ultra-high energies (UHEs) 
by the quantum mechanical decay of supermassive ele- 
mentary "X" particles related to some grand unified the- 
ory (GUT). Such X particles could be released from topo- 
logical defect relics of phase transitions which might have 
been caused by spontaneous breaking of GUT symme- 
tries in the early universe || . TD models of this type are 
attractive because they predict injection spectra which 
are considerably harder than shock acceleration spec- 
tra and, unlike the GZK effect for nucleons, there is no 
threshold effect in the attenuation of UHE 7-rays which 
dominate the predicted flux. 



There has been considerable discussion in the literature 
whether the 7-ray, nucleon, and neutrino fluxes predicted 
by TD scenarios are consistent with observational data 
and constraints at any energy HQ. The absolute flux 
levels predicted by TD models are in general extremely 
uncertain |D3]. Accordingly, in this letter we treat the 
production rate of decaying X particles as a free param- 
eter to be adjusted to match data and constraints. (We 
note that TD scenarios such as annihilation of magnetic 
monopole-antimonopole pairs |fl4| can yield HECR fluxes 
consistent with observations without violating bounds on 
monopole abundances.) Under this "optimal" assump- 
tion, we then use detailed numerical simulations of extra- 
galactic cosmic and 7-ray propagation to show that TD 
models are still viable for an interesting range of param- 
eters. We also explore a signature for TD mechanisms 
based on the isotropic component of the 7-ray flux below 
100 EeV, particularly around 10 EeV. The exposure re- 
quired to test this signature is significantly smaller than 
for measurements above 100 EeV (e.g., as proposed in 
Refs. |1| and |l6|]), as long as discrimination between 7- 
rays and charged cosmic rays (CRs) is possible at a level 
of a few percent. This is within the reach of proposed 
experiments 17 1. 

Top-Down Models. The X particles released, say, in the 
annihilation or collapse of defects such as cosmic strings, 
monopoles, or domain walls could be gauge bosons, Higgs 
bosons, superheavy fermions, etc. depending on the spe- 
cific GUT. Thses X particles would have a mass mx 
comparable to the symmetry breaking scale and would 
rapidly decay typically into a lepton and a quark of 
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roughly comparable energy. We take the primary lep- 
ton produced in a decay to be an electron with energy 
mx /2. (Prior calculations have ignored this lepton which 
is not a good approximation given the lepton's energy.) 
The quark interacts strongly and hadronizes into nucle- 
ons (iVs) and pions, the latter decaying in turn into 7- 
rays, electrons, and neutrinos. Given the X particle pro- 
duction rate, dnx/dt, the effective injection spectrum 
of particle species a (a = 7, N, e , v) via the hadronic 
channel can be written as (dnx /dt)(2/mx)(dN a /dx) > 
where x = 2E/mx, and dN a /dx is the relevant fragmen- 
tation function. For the total hadronic fragmentation 
function dNh/dx we use solutions of the QCD evolution 
equations in modified leading logarithmic approximation 
which provide good fits to accelerator data at LEP ener- 
gies [|18| . We assume that about 3% of the total hadronic 
content consists of nucleons and the rest is produced as 
pions and distributed equally among the three charge 
states. The standard pion decay spectra then give the in- 
jection spectra of 7-rays, electrons, and neutrinos. The X 
particle injection rate is assumed to be spatially uniform 
and in the matter-dominated era can be parametrized as 
dnx/dt oc t~ i+p |J, where p depends on the specific de- 
fect scenario. In this letter we focus on the case p = 1 
which is representative for a network of ordinary cosmic 
string s |19| and annihilation of monopole-antimonopole 
pairs |Lj. 

Numerical Simulations. The 7-rays and electrons pro- 
duced by X particle decay initiate electromagnetic (EM) 
cascades on low energy radiation fields such as the CMB. 
The high energy photons undergo electron-positron pair 
production (PP; 776 — ► e _ e + ), and at energies below 
~ 10 14 eV they interact mainly with the universal in- 
frared and optical (IR/O) background, while above ~ 100 
EeV they interact mainly with the universal radio back- 
ground (URB). In the Klein-Nishina regime, where the 
center of mass energy is large compared to the electron 
mass, one of the outgoing particles usually carries most of 
the initial energy. This "leading" electron (positron) in 
turn can transfer almost all of its energy to a background 
photon via inverse Compton scattering (ICS; e^/b —> e'7). 
EM cascades are driven by this cycle of PP and ICS. The 
energy degradation of the "leading" particle in this cy- 
cle is slow, whereas the total number of particles grows 
exponentially with time. This makes a standard Monte 
Carlo treatment difficult. We have therefore used an an 
implicit numerical scheme to solve the relevant kinetic 
equations. A detailed account of our transport equa- 
tion approach is in Ref. |2^]. We include all EM inter- 
actions that influence the 7-ray spectrum in the energy 
range 10 8 eV < E < 10 25 eV, namely PP, ICS, triplet 
pair production (TPP; ejb — ► ee~e + ), and double pair 
production (77^ — > e~e + e~e + ). The relevant nucleon in- 
teractions implemented are pair production by protons 
(P7& - * pe~e + ), photoproduction of single or multiple 
pions (Njb ~> N nit, n > 1), and neutron decay. Pro- 



duction of secondary 7-rays, electrons, and neutrinos by 
pion decay is also included. We assume a flat universe 
with no cosmological constant, and a Hubble constant 
of Ho = 75 km sec _1 Mpc _1 throughout. An important 
difference with respect to past work is that we follow 
all produced particles, whereas the often-used continu- 
ous energy loss (CEL) approximation (e.g., [pHI ) follows 
only the leading cascade particles. We find that the CEL 
approximation can significantly underestimate the cas- 
cade flux at lower energies. 

Similar studies using somewhat different numerical 
techniques have been performed for the case of dis- 
crete sources injecting 7-rays and nucleons monoenerget- 
ically JTT[] and more recently for fragmentation functions 
oc x~ 1,5 (l — x) 2 and spatially uniform injection ]l2[| . 
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FIG. 1. Predictions for the differential fluxes of 7-rays 
(solid line) and nucleons (long dashed line) by a TD model 
characterized by p = 1, m x = 2 x 10 16 GeV, for no EGMF. 
The dashed line shows the 7-ray flux predicted by the CEL 
approximation. The dash-dotted line shows the result when 
there is no IR/O background. Also shown are the combined 
data from the Fly's Eye § and the AGASA @ experiments 
above 10 EeV (dots with error bars) , piecewise power law fits 
to the observed charged CR flux (thick solid line) and exper- 
imental upper limits on the 7-ray flux at 1 — 10 GeV from 
EGRET data Q (dotted line on left margin). The arrows 
indicate limits on the 7-ray flux from Ref. p3|. 

Results. Fig. [l] shows the results for the 7-ray and 
nucleon fluxes from a typical TD scenario, assuming no 
EGMF, along with current observational constraints on 
the 7-ray flux. The spectrum was normalized in the 
best possible way to allow for an explanation of the ob- 
served HECR events, assuming their consistency with a 
nucleon or 7-ray primary (although a 7-ray primary is 
somewhat disfavored |2l|). The flux below < 20 EeV 
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is presumably due to conventional acceleration and was 
not fit. The shapes of our spectra are similar to those 
of Ref. However, they normalize their spectra to 

match the observed differential flux at 300 EeV, which 
then leads to an overproduction of the integral flux at 
higher energies. We remark that above 100 EeV, the fits 
shown in Figs, [l and ^ have likelihood significances above 
50% (see Ref. [lq] for details) and are consistent with 
the integral flux above 300 EeV esitmated in Refs. 
Since the energy injected at high redshifts is recycled by 
cascading to lower energies, TD models are significantly 
constrained PJlCj] by current limits on the diffuse 7-ray 
background at 1 - 10 GeV Q . Note that the IR/O back- 
ground strongly depletes the background in the range 
10 11 — 10 14 eV, recycling it to energies below 10 GeV (see 
Fig. [l]) . The predicted background is not very sensitive to 
the specific IR/O background model, however |24jj . Con- 
straints from limits on CMB distortions and light element 
abundances from 4 He-photodisintegration are compara- 
ble to the bound from the directly observed 7-rays [ ^o| . 
The scenario in Fig. [j] obeys all current constraints within 
the normalization ambiguities and is therefore quite vi- 
able. 
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FIG. 2. Same as Fig. [j], but for an EGMF of 10~ 9 G. 

Fig. H shows results for the same TD scenario as in 
Fig. [j], but for a high EGMF ~ 10~ 9 G (near the current 
upper limit p6|). In this case, rapid synchrotron cooling 
of the initial cascade pairs quickly transfers energy out 
of the UHE range. The UHE 7— ray flux then depends 
mainly on the absorption length due to pair production 
and is typically much lower ]l^ , p5| . (Note, though, that 
for mx ^ 10 25 eV, the synchrotron radiation from these 
pairs can be above 100 EeV, and the UHE flux is then not 
as low as one might expect.) To match the HECR flux, 



we must inject more X particles, which leads to a factor 
~ 5 increase in the 7— ray background expected below 
10 GeV. For a high EGMF, the constraints from the flux 
limits below 10 GeV are much more severe and in gen- 
eral, TD scenarios in a high EGMF are only marginally 
allowed. 

The energy loss and absorption lengths for UHE nu- 
cleons and photons are short (< 100 Mpc). Thus, their 
predicted UHE fluxes are independent of cosmological 
evolution. The 7— ray flux below ~ 10 11 eV, however, 
scales as the total X particle energy release integrated 
over all redshifts and increases with decreasing p jl0| . For 
mx = 2 x 10 16 GeV, scenarios with p < 1 are therefore 
ruled out (see Figs. [i]and ||), whereas constant comoving 
injection rates (p = 2) are well within the limits. As the 
EM flux above ~ 10 22 eV is efficiently recycled to lower 
energies, the constraint on p is insensitive to mx for low 
EGMF. This in contrast to earlier CEL-based analytical 
estimates ]9|,[l0| • 
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FIG. 3. Predictions for the differential fluxes of 7-rays 
(solid line) and nucleons above 10 15 eV (long dashed line) by 
a uniform, constant comoving density of shock acceleration 
sources up to a redshift of 4, injecting protons with a spec- 
trum ex E~ 2 - 3 up to 10 22 eV, for a vanishing EGMF. 

We now turn to signatures of TD models at UHEs. 
For low EGMF (e.g., Fig. [l]), the full cascade calculation 
predicts 7-ray fluxes below 100 EeV an order of magni- 
tude higher than those obtained using the CEL approxi- 
mation. Again, this shows the importance of non-leading 
particles in the development of unsaturated EM cascades 
at energies below ~ 10 22 eV. Our numerical simulations 
give a 7/CR flux ratio at 10 EeV of ~ 0.1. The experi- 
mental exposure required to detect a 7-ray flux at that 
level is ~ 4 x 10 19 cm 2 secsr, about a factor 10 smaller 
than the current total experimental exposure. These 
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exposures are well within reach of the proposed Pierre 
Auger Cosmic Ray Observatories [|l7| , which may be able 
to detect a neutral CR component down to a level of 
1% of the total flux. In contrast, if the EGMF exceeds 
~ 1CP 11 G, then UHE cascading is inhibited, resulting in 
a lower UHE 7-ray spectrum. In the 10 -9 G scenario 
of Fig. |, the 7/CR flux ratio at lOEeV is ~ 4 x KT 3 , 
significantly lower than for no EGMF. 

Another not well known factor affecting UHE 7-ray 
propagation is the URB for which we used the spectrum 
suggested in Ref. ^7|. A higher overall URB ampli- 
tude correspondingly reduces the 7-ray flux, but with- 
out significantly changing its spectral shape. Thus, as 
long as 7-rays dominate nucleons in the TD component 
above ~ 100 EeV and the total flux is normalized to the 
HECR events, predictions for the 7/CR flux ratio below 
~ 100 EeV are essentially independent of the URB am- 
plitude and the hadronic fraction at injection. An URB 
cutoff frequency lower than 2 MHz (the value we took) 
affects this ratio in a less trivial way with a tendency to 
smaller 7/CR values. 

Fig. H shows spectra resulting from a typical non-TD 
scenario, in this case a uniform distribution of shock ac- 
celeration sources. Such a scenario gives a UHE CR spec- 
trum with a GZK cutoff and 7-rays are only produced as 
secondaries. Our treatment of multiple pion production 
by nucleons leads to secondary 7-ray fluxes somewhat 
higher than in Refs. The key point is that the 

(isotropic) 7/CR flux ratio is < 10~ 3 at 10 EeV, much 
smaller than predicted by TD models in a small EGMF. 
Ratios as high as 10% can only be reached in the direction 
of powerful nearby acceleration sources. The secondary 
7-ray flux generally decreases still further with decreasing 
maximum injection energy and increasing EGMF ]25| ]. 

In summary, some TD-type scenarios for the HECR 
origin are still unconstrained by current data and bounds 
on 7-ray and UHE CR fluxes. For example, if the mean 
EGMF is < 10 -9 G, spatially uniform annihilation of 
magnetic monopoles and antimonopoles is still a viable 
model for GUT scales up to 10 16 GeV. A solid angle av- 
eraged 7/CR flux ratio ~ 10% at ~ 10 EeV would be 
hard to explain by a conventional acceleration origin of 
HECRs, and assuming the TD picture holds, would place 
an independent upper limit of ~ 10~ n G on the poorly 
known EGMF on scales of a few to tens of Mpc. Ab- 
sence of a high 7/CR flux at this level either rules out 
a TD origin for HECRs or implies an EGMF strength 
> 10~ n G. A test of this signature should be possible 
with currently proposed experiments. TD models also 
predict significant neutrino fluxes. Implications of this 
will be considered in a separate publication p9[] . 
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